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в.  с.  САФРОНОВ 
О  РОСТЕ  ПЛАНЕТ  В  ПРОТОН  Л  АНЕТНОМ  ОБЛАКЕ 

Рассматривается  процесс  роста  протопланетных  тел  в  результате  выпадения  на 
них  частиц  газо-пылевого  облака,  вращавшегося  вокруг  Солнца.  Выводится  формула 
для  скорости  роста  планеты  за  счет  малых  частиц,  имеющих  длину  свободного  про- 
бега, малую  по  сравнению  с  расстоянием  от  Солнца.  Формула  с  точностью  до  постоян- 
ного множителя,  равного  двум,  совпадает  с  формулой  О.  Ю.  Шмидта.  Показывается, 
что  космогонически  приемлемые  сроки  роста  Урана,  Нептуна  и  Плутона  могли  быть 
обеспечены  только  при  очень  малых  скоростях  хаотического  движения  частиц  в  этой 
области  на  начальной  стадии  роста.  Оценивается  температура  планет,  обусловленная 
их  разогреванием  ударами  частиц,  падающих  на  поверхность.  Приближенно  оцени- 
вается первоначальное  распределение  температуры  внутри  Земли.  Найдено,  что  ма- 
ксимальная температура  была  на  глубине  около  1000  км  и  составляла  около  1000°. 

В  результате  широкого  обсуждения  теории  происхождения  Земли  и 
планет  О.  Ю.  Шмидта  в  настоящее  время  установилось  более  или  менее 
единодушное  мнение  о  том,  что  планеты  образовались  из  газо-пылевого 
облака,  вращавшегося  вокруг  Солнца.  Но  в  отношении  способа  их  обра- 
.зования  существует  несколько  точек  зрения.  Согласно  одной  из  них  [1] 
планеты  образовались  в  результате  выпадения  вещества  на  имевшиеся 
в  облаке  планетные  «зародыши».  Согласно  другой  [2]  в  тонком  пылевом 
слое  вследствие  гравитационной  неустойчивости  образовалось  множество 
первичных  сгущений  с  массами  порядка  10^^  г  на  расстоянии  Земли  от 
Солнца.  Эти  сгущения,  сталкиваясь  и  соединяясь  друг  с  другом,  выросли 
в  дальнейшем  до  планетных  размеров.  Согласно  третьей  [3]  в  результате 
гравитационной  неустойчивости  сразу  образовались  сгущения  с  массами 
порядка  масс  больших  планет.  Эта  точка  зрения,  казалось  бы,  лучше  всего 
согласуется  с  фактом  существования  систем  спутников  у  Юпитера,  Сатурна 
п  Урана,  расположенных  в  плоскости  экватора  этих  планет.  Однако  она 
встречает  значительные  трудности,  особенно  для  планет,  близких  к  Солнцу, 
так  как  для  наступления  гравитационной  неустойчивости  нужна  высокая 
плотность  вещества  в  облаке,  чтобы  противодействовать  разрушающему 
действию  приливных  сил  гравитационного  поля  Солнца.  В  связи  с  этим 
делается  предположение  [4],  что  масса  протопланетного  облака  была  не 
меньше  массы  современного  Солнца.  Однако  возможность  диссипации 
легких  элементов  из  устойчивых  сгустков  вещества  отвергается  расчетами 
П.  С.  Шкловского  [5].  С  другой  стороны,- естественно  считать,  что  в  сгу- 
щения, образовавшиеся  в  результате  наступления  гравитационной  не- 
устойчивости, вошла  значительная  часть  массы  всего  облака.  Поэтому  масса 
протопланетного  облака,  по  всей  вероятности,  не  могла  во  много  раз  пре- 
восходить современную  массу  всех  планет.  Но  в  таком  случае  необходи- 
мая для  гравитационной  неустойчивости  высокая  плотность  могла  воз- 
никнуть лишь  в  сильно  уплощенном  диске  толщиной  порядка  размеров  пла- 
неты и  на  два-три  порядка  меньшей,  чем  ширина  зоны  вещества,  вошедшего 
в  состав  планеты.  Из  такого  тонкого  и  широкого  кольца  не  могло  сразу 
образоваться  одно  сгущение. 

Выбор  между  первой  и  второй  точками  зрения  может  быть  сделан  лишь 
после  расчета  скорости  уплощения  пылевой  составляющей  облака  и  ско- 
рости роста  в  нем  планетных  зародышей. 

В  настоящей  статье  исследуются  некоторые  закономерности  роста 
планетных  тел  в  результате  выпадения  на  них  вещества  из  газо-пылевого 
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облака.  Случай  наступления  гравитационной  неустойчивости  здесь  не- 
рассматривается. 

В    статье   не   рассматриваются   процессы  дробления   сталкивающихся 
тел.  Если  отвлечься  от  возможности  частичного  разрушения,  то  это  приблп 
жение   можно   считать   допустимым, —  рассматривается   рост    таких   тел 
которые  вырастают  в  планеты,  не  испытав  ни  одного  столкновения,  приво- 
дящего к  их  разрушению.  По  существу  только  такие  тела  для  нас  и  пред- 
ставляют непосредственный  интерес. 

Для  рассмотрения  процесса  эволюции  протопланетного  облака  нужно 
знать  состояние  облака  в  некоторый  «начальный»  момент.  Однако  в  этом 
вопросе  существует  значительная  неопределенность.  В  частности,  до  сих 
пор  неясно,  образовалось  ли  это  облако  совместно  с  Солнцем,  или  было  им 
захвачено.  В  основу  настоящей  работы  положены  следующие  предположе- 
ния: 

1)  Первоначальная  масса  облака  не  превосходила  существенно  совре- 
менную массу  планет.  Химический  состав  облака  был  аналогичен  совре- 
менному химическому  составу  планет. 

2)  Облако  первоначально  было  достаточно  уплощенным  и  напоминало 
сплюснутый  эллипсоид  вращения  с  отношением  полуосей  порядка  1  :  10 
или  менее.  Вещество  в  любом  элементарном  объеме,  выделенном  в  облаке, 
двигалось  вокруг  Солнца  по  круговой  орбите  и  с  круговой  скоростью  на 
данном  расстоянии  от  Солнца. 

3)  Заметная  доля  вещества  облака  находилась  в  пылевом  состоянии 
(во  всяком  случае  не  меньшая,  чем  в  темных  туманностях).  Поэтому  облако 
было  непрозрачным  для  солнечного  излучения  и  относительно  холодным. 
Атомы,  сталкивавшиеся  с  пылинками,  конденсировались  на  них  после 
небольшого  числа  столкновений.  В  результате  большинство  атомов  должно 
было  быстро  осесть  на  пылинках.  Следовательно,  ведущая  роль  в  процессе 
эволюции  облака  принадлежала  пылевым   частицам. 

4)  В  облаке  имелись  тела  («зародыши  планет»),  достаточно  массивные 
для  того,  чтобы  пылевые  частицы,  сталкивавшиеся  с  ними  с  малыми  отно- 
сительными скоростями,  удерживались  на  поверхности   этих  тел. 

Полученные  ниже  результаты  справедливы  лишь  в  рамках  этих  пред 
положений. 

Формула  для  скорости  роста    планет 

В  1945  г.  О.  Ю.  Шмидт  [1]  предложил  следующую  формулу  для  скорости 
роста    планет 

Где  Яу  ж  Нч  —  расстояния  от  Солнца  внутренней  и  внешней  границ  «зо- 
ны питания»  данной  планеты,  т.е.  области,  в  которой  двигались  частицы  ,^ 
вошедшие  в  дальнейшем  в  состав  планеты;  ^  —  первоначальная  масса 
вещества  в  этой  зоне;  Р  —  период  обращения  планеты  вокруг  Солнца; 
т  ж  г  —  масса  и  радиус  планеты  в  момент  I. 

При  выводе  этой  формулы  предполагалось,  что  частицы  движутся  по 
кеплеровым  орбитам  вокруг  Соднца.  Два  раза  за  период  своего  обраще- 
ния вокруг  Солнца  они  пересекают  плоскость  орбиты  растущей  планеты, 
и  часть  из  них  при  этом  сталкивается  с  планетой  и  остается  на  ней.  Грави- 
тационное воздействие  планеты  на  частицу  не  учитывалось. 

Если  современную  массу  планет  рассеять  по  сравнительно  тонкому 
плоскому  слою  (толщиной  в  0.01—0.1  от  его  диаметра),  простирающемуся 
до  границ  современной  солнечной  системы,  то  плотность  в  этом  слое  ока- 
зывается порядка  10"^^  г/см^.  Это  значение  принимается  в  ряде  работ.  При 
такой  плотности  вещества  в  протопланетном  облаке  по  кеплеровым  орбитам 
могут  двигаться  только  наиболее  крупные  частицы  с  радиусом  в  несколько 
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сантиметров  и  более.  Частицы  меньших  размеров  (  и  тем  более  атомы  и  мо- 
лекулы газа),  участвуя  в  общем  движении  вокруг  Солнца,  кроме  того, 
часто  сталкиваются  друг  с  другом.  Их  движения  отличны  от  движения  по 
кеплеровым  орбитам.  Повидимому,  первоначально  протопланетное  облако 
состояло  в  основном  из  мелких  частиц.  Следовательно,  необходимо  опре- 
делить скорость  роста  планеты  в  результате  выпадения  мелких  частиц, 
хаотически  движущихся    относительно    нее. 

Выражение  для   скорости  роста  можно  записать  в    следующем  виде: 

Ж  =  -^Г-  (2) 

Здесь  р  —  плотность  вещества  в  облаке,  V  —  средняя  (по  плотности) 
скорость  частиц  относительно  рассматриваемого  тела,  т  —  масса  тела. 
Для  не  слишком  больших  тел,  когда  их  гравитационным  действием  на 
частицы  можно  пренебречь,  г  представляет  геометрический  радиус  тела. 
Для  больших  тел  (начиная  с  размеров  крупных  астероидов)  под  г  следует 
понимать   эффективный  радиус. 

Формула  (2)  использовалась  многими  авторами.  Вейцзеккер  [6],  Хойл 
[7],  Тер-Хаар  [8]  применяли  ее  для  подсчета  скорости  конденсации  газа 
на  планетном  зародыше.  В  эту  формулу  входят  неизвестные  величины  р  и 
V,  которые  исключены  из  формулы  О.  Ю  Шмидта.  Плотность  р  по  мере 
роста  планеты  уменьшается.  Весьма  неопределенна  величина  скорости  час- 
тиц V,  так  как  неизвестно,  какие  скорости  преобладали  в  облаке  —  тер- 
мические или  турбулентные.  Принятие  в  работах  [6 — 8]  некоторых  наугад 
взятых  значений  скорости  V  делало  оценки  скорости  роста  весьма  грубыми. 
Принятие  же  постоянного  значения  плотности  р  приводило  к  неправиль- 
ному выводу  о  высокой  первоначальной  температуре  поверхностных  слоев 
Земли. 

Так  как  и  вещество  в  облаке,  и  планеты  предполагаются  движущимися 
вокруг  Солнца  с  круговой  скоростью,  средняя  скорость  частиц  V  относи- 
тельно планеты  равна  их  хаотической  скорости.  Хаотическая  скорость  не- 
посредственно связана  с  толщиной  протопланетного  диска,  а  через  нее  — 
€  плотностью  вещества  в  облаке.  Поэтому  можно  попытаться  выразить  V 
■через  р  и  упростить  выражение  (2). 

Если  обозначить  через  Н  толщину  однородного  кольца,  имеющего  мас- 
су ^  —  ти  плотность,  равную  максимальной  плотности  ро  в  центральной 
плоскости,  в  которой  расположена  орбита  планеты,  то,  очевидно,  имеет 
место   соотношение: 

ртг  {Н1  ~Я\)Н  =  д  —  т.  (3) 

Выражение  (2)  можно  теперь  переписать  в  виде: 

Лт  г^ 


сИ         Г  7?2  _  р2 


{В'^-Я^)Н 


{^-т)V.  (4) 


Для  определения  зависимости  между  V  ш  Н  применим  к  частицам  соот- 
ношения кинетической  теории  газов.  Условие  гидростатического  равнове- 
сия можно  записать  в  виде  «» 

^^  =  -Р2-  (5) 

Давление  р  равно 

^^  =  ТР^^'  (6) 

1де  V^  —  средний  квадрат  скорости  частиц  в  облаке.  Составляющая  силы 
тяготения  по  2-координате  ^  складывается  из  двух — силы  тяготения, 
вызванной  Солнцем,  и  силы  тяготения  самого  диска.    Вторая  становится 
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сравнимой  с  первой  только  в  случае  очень  сильно  уплощенного  облака, 
когда  в  нем  возникает  гравитационная  неустойчивость.  В  рассматриваемом 
нами  случае  менее  плоского  облака  ею  вполне  можно  пренебречь.  Тогда 

ат.  лз     •  \ч 

Предположим,  что  все  частицы  одного  размера  и  что  ч)  не  зависит 
от  высоты  2,  т.  е.  что  кинетическая  температура  одинакова  по  всей 
толщине  диска.  Тогда  уравнение  (7)  легко  интегрируется  и  дает 

ЗСМ0   ,, 

р(2)  =  рое"^^^    .  (8) 

Следовательно,  толщина  однородного  диска  равна 

— оо 

Подставим  это  значение  Н  в  формулу  (4)  и  заменим  V^  на  1^^,  .  со- 
гласно известному  соотношению  для  максвелловского  распределения  ско- 
ростей   V^  =  -о-^'^.    Мы   получим 


ат  2  ,  /<^^^©  г^  ,^  . 


(10) 


Так  как  период  обращения  планеты  вокруг  Солнца  равен  Р  = 
=  27Г  I/  -^пгг ,  то  формулу  (10)  можно  записать  в  виде 

Таким  образом  найденная  нами  формула  для  скорости  роста  планет  с 
точ^^остью  до  постоянного  множителя,  равного  двум,  совпадает  с  форму- 
лой О.  Ю.  Шмидта   (1). 

На  первой  стадии  роста  тела,  когда  его  масса  еще  мала,  влиянием 
притяжения  этого  тела  на  движение  сталкивающихся  с  ним  частиц 
можно  пренебречь.  В  этом  случае  г  в  формулах  (1)  и  (И)  представляет 
просто  геометрический  радиус  тела,  предполагаемого  сферическим.  При 
большой  массе  тело  заметно  притягивает  частицы,  что  увеличивает  ча- 
стоту их  столкновений  с  телом.  Расчет  различен  для  крупных 
частиц,  не  сталкивающихся  друг  с  другом,  и  для  мелких  частпц. 
или  атомов  и  молекул  газа.  В  первом  случае  рассматривается  движение 
одной  изолированной  частицы  в  поле  тяготения  тела.  Вследствие  искри- 
вления траектории  частицы  в  направлении  притягивающего  тела  эффектив- 
ное сечение  ее  столкновения  с  этим  телом  увеличивается: 

и/2  =  ^.2^1  +  ^).  (12) 

Здесь  т  ж  г  —  масса  и  радиус  тела,  V  —  скорость  частицы  относительно 
тела  «в  бесконечности».  Будем  для  краткости  называть  величину  I  радиу- 
сом столкновения.  Во  втором  случае  необходимо  уже  газодинамическое 
рассмотрение  движения  сплошной  среды  в  поле  тяготения  сферического 
тела.  Точное  решение  этой  задачи  никем  не  производилось.  Рассмотрен 
частный  случай  [9],  когда  тело  движется  относительно  среды  с  некоторой 
скоростью  V  и  когда  в  этой  среде  имеется  достаточная  примесь  частиц,  ко- 
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торые  уменьшают  скорости  беспорядочного  движения  частиц  среды  вслед- 
ствие неупругих  столкновений  и  делают  возможным  захват  частиц  телом. 
Радиус  захвата  г*  при  этом  процессе,  названном  аккрецией,  выражается 
в    виде 

'•  =^;г-.  (13) 

граница  области  захвата  характеризуется  равенством  средней  кинетиче- 
ской энергии  движения  частиц  относительно  тела  и  потенциальной  энер- 
гии частиц  в  ноле  тяготения  этого  тела.  По  существу  формула  (13)  дает 
лишь  верхний  предел  для  радиуса  захвата.  Повидимому,  действительное 
значение  г*  должно  быть  меньше  и  определяется  условиями  конкретной 
задачи  (плотностью  газа  и  частиц,  степенью  неупругости  столкновений, 
скоростью  движения  среды  относительно  тела). 

Формула  (13)  была  затем  обобщена  на  случай,  когда  тело  неподвижно 
относительно  среды.  При  этом  под  V  стало  подразумеваться  значение  хао- 
тической скорости  частиц  или,  что  почти  то  же,  скорости  звука. 

Имеется  еще  третий,  промежуточный  случай,  когда  плотность  частиц 
еще  недостаточна  для  того,  чтобы  необходимо  было  газодинамическое 
рассмотрение  движения  сплошной  среды,  но  уже  достаточна  для  того,  чтобы 
частицы  испытывали  многократные  столкновения  друг  с  другом  в  области 
повышенной  плотности  на  линии  движения  тела  относительно  облака. 
Подобный  случай  был  исследован  Т.  А.  Агекяном  [10]  при  рассмотрении 
возможности  захвата  вещества  звездой,  проходящей  через  пылевое  облако. 

Легко  найти  значение  т,,,  при  котором  радиус  захвата  становится 
равным  геометрическому  радиусу 

где  8  —  средняя  плотность  тела. 

В  выражении  (12)  при  этом  значении  массы  второе  слагаемое  в  скобках 
в  правой  части  становится  равным  единице.  Можно  принять,  что  при  массе 
тела,  меньшей  т^,  притяжением  тела  можно  пренебречь  и  в  (11)  использо- 
вать обыкновенный  геометрический  радиус.  При  гпГ^то  для  крупных  ча- 
ст'иц  следует  брать  радиус  столкновения  /,  а  для  мелких  частиц  и  газа — 
радиус  захвата  г*.  Из  (14)  видно,  что  величина  то  существенно  зависит 
от  скорости  частиц.  При  8  =  3  и  г;  =  1  км/сек  получим  то  =  5.8'  10^^  ги  го  = 
=  780  км. 

Время  роста  и  температура  планет 

В  соответствии  с  изложенными  в  начале  статьи  соображениями  о  ха- 
рактере эволюции  протопланетного  облака  и  образования  в  нем  планет  мож- 
но произвести  приблизительную  оценку  продолжительности  отдельных 
этапов  роста  планет. 

Начальная  стадия  роста.  На  первой  стадии  раста  тела,  когда 
его  масса  мала,  и  потому  мало  его  гравитационное  воздействие  на 
частицы,  радиус  захвата  при  аккреции  г*  оказывается  меньше  геомет- 
рического радиуса  тела.  Однако  фактический  радиус  захвата  при  этом 
с  достаточно  хорошей  точностью  определяется  выражением  (12),  незави- 
симо от  размеров  частиц,  сталкивающихся  с  телом.  Подставив  в  (И) 
вместо  г  значение  радиуса  столкновения  /  согласно  (12),  получим 


Это  выражение  легко  привести  к  виду 

с1г  =  А{1  —  г^)  (1  +  1г^)  сИ,  (16) 
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где 


Интегрируя  (15),  получим 


(17) 


1 ).  Ш  [1  +  1{^у^]  -  ^_  (1  _  Х^)  +  ^(1  -  л^)  агс  Щ  - 


Н-Г/.агс1^уТ(-^)'^^|. 


+ 
(18) 


При  X  >  1  и  при  значениях  т,  не  слишком  близких  к  ^,  в  (18)  прео- 
бладающим является  последний  член.  В  этом  случае  с  достаточно  хоро- 
шей точностью  можно  написать 


А1  = 


1 


^^^Vс^§V}:^^^ 


т  у/з 


(19) 


Это  же  выражение  получается  при  интегрировании  (15),  если  в  нем  вме- 
сто ^ — т  написать  ^.  Таким  образом  оно  годится  для  начальной  стадии 
роста  планеты,  когда  т  <  ^.  Отсюда  можно  сразу  найти  значение  1{то). 
Так  как 


то 


Кг(^)-  =  1, 


г  {шо)  = 


4ЛГХ 


(20) 
(21) 


Приняв  некоторые  значения  Я^  и  Лг  для  границ  зон  планет  так,  чтобы 
плотность  в  облаке  менялась  от  планеты  к  планете  более  или  менее  плавно, 
и  взяв,  например,  8  =  3  для  внутренних  планети8  =  1.3  для  внешних  планет, 
можно  найти  приблизительные  значения  Л  и  X  для  каждой  планеты,  а  затем 
время  роста  до  массы  то.  Найденные  таким  образом  величины  -4,  X  и  /  (то) 
приведены  в  табл.  1  для  скорости  частиц  V  =  1  км/сек.  Масса  Плутона 
условно  принята  равной  массе  Земли.  Легко  найти  время  роста  для  дру- 
гих значений  скорости  V,  так  как  это  время  пропорционально  V. 

Таблица    1 


Планета 


А,  сек.-' 


/  ('Л,),  лет 


Меркурий 
Венера  . 
Земля     . 
Марс  .    .    , 
Юпитер 
Сатурн  . 
Уран  .   .    , 
Нептун  . 
Плутон 


10-16 

10-16 
10-16 
10-17 

10-" 
10-18 
10-20 
10-20 

10-21 


15 

1 

90 

9  • 

102 

2  • 

23 

5  • 

3600 

3  . 

1610 

3  • 

460 

3  . 

520 

1  . 

78 

1  • 

10^ 
106 
10' 
108 
10' 
108 
1010 

1011 
1012 


Из  рассмотрения  этой  таблицы  видно,  что  рост  внешних  планет,  начи- 
ная с  Урана, не  укладывается  в  космогонически  приемлемые  сроки  при  тех 
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скоростях  V,  которые  принимаются  как  наиболее  вероятные  в  работах 
[6 — 8].  Процесс  роста  мог  протекать  быстрее  лишь  при  меньших  значениях 
скорости  частиц  в  облаке.  Для  того  чтобы  продолжительность  стадии  роста 
планеты  от  нулевой  массы  до  массы  то  была  равной,  например,  Ю**  лет, 
скорость  частиц  должна  быть  равной  в  зоне  Урана  10^  см/сек,  в  зоне  Не- 
птуна —  4  •  10^  см/сек,  в  зоне  Плутона  всего  лишь  20  см/сек.  Для  того  чтобы 
п  зоне  этих  планет  могла  наступить  гравитационная  неустойчивость,  необ- 
одимы  скорости  частиц,  на  порядок  меньшие. 
Отсюда  следуют  два  вывода:  а)  планетные  зародыши  в  зоне  Урана — 
Плутона  не  могли  вырасти  до  массы  то  в  газовом  облаке  в  результате 
конденсации  на  них  атомов  газа,  как  это  предполагается  в  работах  Хойля 
и  Вейцзеккера.  Только  выпадение  пылевых  частиц,  двигавшихся  с  малыми 
скоростями,  могло  обеспечить  эффективный  рост  этих  зародышей;  б)  тур-^ 
булентные  движения  в  зоне  внешних  планет  не  играли  на  первой  стадии 
существенной  роли,  так  как  в  противном  случае  частицы  имели  бы  значи- 
тельные скорости,  и  сроки  роста  планет  до  масы  то  были  бы  неприемлемо 
велики. 

Из  формулы  (И)  видно,  что  скорость  роста  планет  на  первой  стадии  не 
зависит  от  скорости  частиц  V  и  может  быть  вычислена  без  дополнительных 
предположений.  Можно,  например,  найти  время  роста  тел  до  размеров 
первичных  сгуш;ений,  если  их  массы  меньше  то.  В  табл.  2  приведены  массы 
первичных  сгущений  т^,  найденные  по  формуле  Л.  Э.  Гуревичаи  А.  И.  Ле- 
бединского [2],  и  время  роста  планет  до  этой  массы. 


Таблица    2 


Планеты 

Меркурий 

Венера 

Земля 

Марс 

т^   (грахМм) 

1  (т^)   (лет)     .... 

10» 
7  .  102 

4  .  101* 
4.  103 

1015 
2  .  10* 

1012 

105 

Из  рассмотрения  таблицы  видно,  что  для  внутренних  плант  масса 
первичных  сгущений  достигалась  в  результате  роста  тел  очень  быстро. 
Это  обстоятельство  наряду  с  наличием  в  этой  области  больших  прилив- 
ных сил,  а  также  более  высокой  температуры,  затрудняющих  возникно- 
вение условий  для  гравитационной  неустойчивости,  свидетельствует  в  поль- 
зу механизма  роста  внутренних  планет  (Меркурия,  Венеры,  вероятно. 
Земли  и,  может  быть.  Марса)  путем  выпадения  вещества  на  первоначально 
совсем  малые  планетные  зародьппи. 

Рост  болеекрупных  тел.  После  того  как  масса  тела  достигла  зна- 
чения то,  в  формуле  (11)  в  качестве  г  необходимо  уже  брать  не  геометри- 
ческий радиус,  а  либо  радиус  столкновения  I,  либо  радиус  захвата  г*.  Рас- 
смотрим оба  эти  случая  отдельно. 

а)  Скорость  роста  тела  в  результате  выпадения  на  него  крупных 
частиц,  не  сталкивающихся  друге  другом,  определяется  соотношением  (15). 
Время  роста  от  нулевой  массы  до  массы  т  определяется  формулой  (18)  или 
приближенной  формулой  (19).  При  т,  приближающихся  к  ^,  из  (19)  и  (21) 
находим 


Ц^)^2({т,). 


(22) 


Таким  образом  продолжительность  роста  тела  до  массы  то  приблизительно 
равна  продолжительности  роста  от  то  до  т,  близкого  к  ^,  например,  до 
т  =  0.9  ^  или  до  т  =  0.99  ^.  Однако  в  связи  с  исчерпанием  вещества  в 
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облаке  скорость  дальнейшего  процесса  роста  асимптотически  замедляется 
до  нуля.  Поэтому  мы  не  можем  взять  т  =  ^,  так  как  в  этом  случае  выра- 
жение (18)  расходится  и  переход  к  (19)  незаконен. 

Если,  например,  мы  примем  V  =  {  км/сек,  то  для  Земли  л  =  10^.  В  этом 
случае  при  т  =  ^  (1 — 10~^)  логарифмические  члены  в  (18)  равны  члену,  со- 
храненному в  (19),  следовательно,  время  роста  до  этой  массы  в  два  раза  боль- 
ше, чем  получается  из  (19).  Поскольку  полученные  значения  для  времени 
роста  на  второй  стадии  практически  повторяют  значения  табл.  1,  то,  как  и 
прежде,  чтобы  получить  космогонически  приемлемые  сроки  роста  внешних 
планет,  необходимо  принять,  что  скорости  частиц  в  облаке  были  значи- 
тельно меньше  обычно  принимаемых  скоростей,  порядка  1  км/сек.  Малая 
скорость  частиц  необходима  в  начале  второй  стадии,  когда  процесс  роста 
из-за  малой  массы  планеты  развивается  медленно. 

б)  Количество  выпадающих  на  тело  мелких  частиц,  длина  свободного 
пробега  которых  не  превышает  значительно  размеров  тела,  определяется 
радиусом  захвата.  В  первом  приближении  радиус  захвата  можно  выра- 
зить в  виде  (13).  При  этом  предполагается,  что  столкновения  между  части- 
цами (или  частицами  и  атомами)  обладают  достаточной  степенью  неупру- 
гости, приводящей  к  уменьшению  их  хаотических  скоростей  в  области  по- 
вышенной плотности  частиц  на  линии  движения  тела  относительно  облака 
и  вследствие  этого  к  возможности  захвата  частиц. 

Подставив  в  (И)  значение  г*  согласно  (13),  получим 


^  =  Вт^{^-т),  (23) 


где 


о               16С2  . 

В  =  9 9 •  (24) 

Интегрируя  (23),  находим  время  роста  планеты  от  массы  т^  до  массы  т 


В1 


=^[1.+'Ч^.-^)-1-Ч-^-0]-         <^^) 


Величина  то,  как  и  ранее,  определяется  из  (14).  В  пра;вой  части  (25) 
преобладающим  является  первый  член.  Так  как  т^^  —  V^,  а  ^-^'у"*,  то 
время  роста  приблизительно  пропорционально  скорости  частиц  1)  отно- 
сительно планеты.  При  т^  <^  ^  рост  планеты  путем  аккреции  от  массы 
то  до  современной  массы  происходит  в  3 — 4  раза  быстрее,  чем  ее  рост 
в  начальной  стадии  от  нулевой  массы  до  т^. 

Первоначальная  температура  планет.  При  ударе  частиц 
о  поверхность  планеты  их  кинетическая  энергия  превращается  в  тепло, 
которое  частично  идет  на  нагревание  поверхностного  слоя,  частично 
излучается  в  пространство.  С  достаточной  для  наших  расчетов  степенью 
точности  температуру  поверхности  мвжно  определить  следующим  соот- 
ношением [И]: 

^!1^  =  47гг2аГ4.  (26) 

Здесь  написаны  главные  члены  в  общем  соотношении  баланса  энергии: 
слева  —  потенциальная  энергия,  освобождающаяся  при  выпадении  вещества 
на  поверхность  планеты;  справа  —  энергия,  излучаемая  поверхностью 
планеты.  Члены,  характеризующие  энергию,  идущую  на  нагревание 
поверхности,  на  фазовое  превращение  падающих  частиц,  кинетическую 
энергию  частиц  до  их  сближения  с  планетой,  малы  по  сравнению  с  на- 
писанными выше  и  потому  исключены. 


м 
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50? 


Выражение  (26)  можна  переписать  в  виде 


гр,  ио  ат 


ОВОт 
Зст   (II 


(27) 


Температура  поверхности  растущей  планеты   определяется  величиной 

с1т  с.  Т-,  ■       ^ 

-7- ,  Т.  С.  интенсивностью  процесса  роста  в  данный  момент.  В  начальной 


с11 


(1т 


стадии    роста    -^~  очень    мало,  и   потому  температура  тоже  мала.    При 

-,  „  „     (1т 

массе -планеты,  близкой  к   современной,  -у-  тоже   мало,    вследствие  ис- 
черпания   вещества    в  облике    [мал    множитель    {^  —  т)    в  (И)].  Легко 

найти  значение  т,  при   котором  -г-   достигает    максимума.    Для    выра- 

4 
жения  (15)  этот  максимум  будет   при    т  =  -^  (^,     для     выражения     (23) 

максимум    достигается   при  т  =  -^^. 

Соответствующие    этим  величинам  т    максимальные  значения  темпе- 
ратуры для  обоих  случаев  приведены  в  табл.  3. 

Таблица    3 


Планеты 

Мер- 
курий 

Вене- 
ра 

Земля 

Марс 

Юпи- 
тер 

Сатурн 

Уран 

Неп- 
тун 

Плу- 
тон 

т(■^^) 

г(4<?) 

190 
420 

630 
1700 

510 
1400 

80 
150 

'2700 
17000 

660 
3600 

140 
600 

120 
520 

25 
75 

Расчет  проведен  при  скорости  частиц  V  =  1  км/сек  и  при  значениях  о, 
несколько  меньших  средней  плотности  планет.  Первая  строка  рассчитана 
для  случая,  когда  на  планету  падают  крупные  частицы,  вторая  —  когда 
падают  мелкие.  Так  как  в  протопланетном  облаке  существовали,  невиди- 
мому, и  те  и  другие,  то  действительное  значение  Т  должно  быть  заключено 
где-то  между  приведенными  выше.  Результат  существенно  зависит 
от  принятого  значения  г;.  В  первом  случае  Т-^г;— '/%  во  втором  —  Т'-^у^. 
Уверенные  значения  для  температуры  планет  можно  будет  получить  лишь 
после  уточнения  величины  скорости  частиц  в  облаке. 

Как  уже  указывалось  вьппе,  значение  V  =  1  км/сек  неприемлемо  для 
начальной  стадии  роста  Урана,  Нептуна  и  Плутона.  Однако,  если  значения 
V,  приемлемые  для  начальной  стадии,  подставить  в  (27),  то  в  случае  ак- 
креции температура  Урана,  Нептуна  и  Плутона  на  последующей  стадии 
роста  получается  слишком  высокой  —  несколько  десятков  тысяч  градусов. 
Л.  Э.  Гуревичем  и  А.  И.  Лебединским  [2]  было  отмечено,  что  первичные 
сгущения  в  результате  гравитационного  взаимодействия  должны  были 
приобретать  значительные  относительные  скорости.  Из  полученных  выше 
результатов  следует,  что  приемлемые  значения  для  времени  роста  и  темпе- 
ратуры планет  можно  получить  лишь  при  условии,  если  скорости  частиц 
относительно  планеты  на  второй  стадии  роста  были  значительно  больше, 
чем  на  первой.  Таким  образом  из  совершенно  различных  и  независимых 
соображений  следует  один  и  тот  же  вывод  об  изменении  относительных  ско- 
ростей частиц  в  облаке  в  процессе  роста  планет.  При  отсутствии  в  облаке 
крупных  тел  скорости  частиц  уменьшались,  и  облако  уплощалось.  По  мере 
образования  крупных  тел  скорости  в  облаке  стали  расти,  а  облако  утол- 
щаться. 
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Первоначальное  распределение  температуры  внут- 
ри Земли.  Условимся  понимать  под  выражением  «первоначальный» 
такой  момент  времени,  когда  Земля  уже  практически  достигла  совре- 
менных размеров,  но  момент,  настолько  удаленный  в  прошлое,  что  время 
роста  Земли  до  этого  момента  времени  можно  считать  малым  по  сравнению 
со  временам,  протекшим  после  этого  момента.  Ввиду  медленности  про- 
цесса разогрева  Земли  при  распаде  радиоактивных  элементов  такое  опре- 
деление позволяет  пренебречь  влиянием  радиоактивного  распада  на  терми 
ческий  режим  Земли  на  первом  этапе  ее  формирования. 

На  основании  приведенных  выше  данных  нетрудно  найти  перво- 
начальное, распределение  температуры  внутри  Земли  по  радиусу. 

Результаты  изображены  на  рис.  1  соответственно  для  случаев  выпадения 
на  Землю  крупных  тел  и  мелких  частиц.  По  оси  абсцисс  отложена  вели- 
чина т'/»  в  единицах  современной  массы  Земли,  которая,  очевидно, 
совпадает  с  г  в  случае  постоянной  плотности.  Расчет  проведен  для  ско- 
ростей частиц  V  =  0,5;  1  и  3  км/сек.  Средняя  плотность  Земли  принята 
меняющейся  от  3  при  т  =  0.01  (?  до  5.52  при  т  =  <?•  Если  все  приобретен- 


е  0,1  0,2  0,0 0,^  0,5  0,00,7  0,8  0,9  1,0  т''' 


О   О,}  0,2 0,3 0,^  0,5  0,6 0,7  0,8 0,9 [О  т^'^ 


Рис.  1.  Первоначальное  распределение  температуры  внутри  Земли.  Слева — рост 
Земли  происходил  в  результате  выпадения  на  нее  крупных  тел,  справа — в  ре- 
зультате выпадения  мелких  частиц  (аккреция) 


ное  таким  образом  количество  тепла  равномерно  распределить  по  всей 
Земле,  то  получаются  соответственно  следующие  средние  температуры:  в 
первом  случае  Т  =630,  460  и  ЗЭ0°К,  во  втором  случае  Т  =  2400,  1200  и 
450°  К.  Можно  полагать,  что  уже  на  ранней  стадии  роста  Земли  облако 
стало  прозрачным  и  частицы  в  облаке  приобрели  температуру  совре- 
менных окресностей  Земли   (-^280°). 

Здесь  не  учтено  нагревание  Земли  при  ее  сжатии  под  давлением 
прибавляющихся  сверху  новых  слоев.  По  некоторым  данным  эта  величина 
составляет  несколько  сот  градусов   [12]. 

Увеличение  в  процессе  роста  планеты  скорости  частиц  означает,  что 
реальному  первоначальному  распределению  температуры  внутри  Земли  не 
соответствует  на  рис.  1  ни  одна  из  кривых,  вычерченная  для  какого- 
либо  постоянного  значения  V.  Точки  кривой  фактического  распределения 
должны  переходить  постепенно  с  одной  кривой  на  другую  от  малых  зна- 
чений скорости  г;  на  первой  стадии  к  большим  — на  второй  стадии.  Только 
знание  действительного  значения  V  как  функции  времени  позволит  полу- 
чить истинное  распределение  первоначальной  температуры  внутри  Земли, 

Так  как  найденные  значения  температуры  существенно  зависят  от  пр: 
пятых  значений  V,  то  желательно  было  бы  оценить  V  хотя  бы  при  близ  и*^ 
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тольно.  Выше  предполагалось,  что  вещество  в  облаке  не  имеет  системати- 
ческой скорости  относительно  планеты.  Однако,  как  известно,  Земля  дви- 
жется не  по  круговой  орбите,  а  по  эллиптической  с  эксцентриситетом 
с  =  0.017.  Поэтому  относительно  вещества,  движущегося  по  круговой  орбите 
с  тем  же  периодом  обращения,  Земля  имеет  систематическую  скорость 
порядка  Уорб  б,  т.  е.  0,5  км/сек.  Это  отклонение  от  круговой  орбиты  является 
результатом  возмущающего  действия  на  Землю  со  стороны  других  планет. 
\[ожно  думать,  что  в  процессе  роста  Земля  подвергалась  возмущениям  та- 
кого же  порядка  как  со  стороны  других  планет,  так  и  со  стороны  тел,  во- 
шедших затем  в  состав  Земли.  Тела  же,  падавшие  на  Землю,  испытывали, 
повидимому,  еще  большие  возмущения  в  гравитационном  поле  Земли. 
Поэтому  можно  полагать,  что  скорости  тел  и  частиц  относительно  Земли 
были  не  менее  1  км/сек.  Из  рассмотрения  рис.  1  можно  сделать  заключе- 
ние, что  температура  внутри  Земли  первоначально,  повидимому,  не  прево- 
сходила 1000°  К  и  была  максимальной  на  глубине  порядка  1000 — 1500  км. 

Температура  самых  поверхностных  слоев  Земли,  после  того  как  ее 
рост  в  основном  закончился,  никогда  не  была  высокой  ввиду  исключитель- 
ной медленности  процесса  роста  на  заключительной  стадии. 

Таким  образом  вытекающее  из  теории  О.  Ю.  Шмидта  представление  о 
первоначальной  температуре  Земли  совершенно  не  похон^е  на  существовав- 
шее ранее  представление  об  огненно-жидкой  Земле,  отделившейся  от  Солн- 
ца и  имевшей  температуру,  близкую  к  температуре  внешних  слоев  Солнца, 
т.  е.  порядка  5000° К. 

В  этом  именно  смысле  и  следует  понимать  употребляемое  О.  Ю.  Шмид- 
том выражение,  что  Земля,  образовавшаяся  из  газо-пылевого  облака,  пер- 
воначально была  холодной.  Температура  внутри  Земли  могла  достигать 
тысячи  градусов.  Однако  со  времени,  когда  рост  Земли  в  основном  закон- 
чился (например,  начиная  с  т  =  0,99  0,  поверхность  Земли  всегда  была 
действительно  холодной. 

Заключение 

Рассмотрение  процесса  роста  протопланетных  тел  в  результате  выпаде- 
ния на  них  частиц  показывает,  что  при  помощи  этого  механизма  можно 
объяснить  образование  планет  из  газо-пылевого  облака  в  космогонически 
приемлемые  сроки.  Рассмотрение  процесса  роста  внешних  планет  —  Урана, 
Нептуна,  Плутона  —  приводит  к  выводу,  что  скорости  частиц  в  облаке 
первоначально  были  очень  малыми,  а  затем  по  мере  их  роста  увеличива- 
лись. К  сожалению,  незнание  действительных  значений  этой  скорости  при- 
водит к  значительной  неопределенности  численных  оценок  скорости  роста 
планет  и  первоначальной  температуры  планет.  В  отношении  Земли  можно 
лишь  сказать,  что  температура  слоев,  расположенных  на  глубине 
1000 — 1500  км,  могла  быть  порядка  тысячи  градусов.  Однако  самые  поверх- 
ностные слои  Земли  (уже  в  основном  достигшей  современных  размеров) 
всегда  были  сравнительно  холодными. 

Для  получения  надежных  количественных  оценок  необходимо  решение 
ряда  вопросов,  касающихся  эволюции  протопланетного  облака.  Совершен- 
но необходимо  исследование  скорости  диссипации  различных  химических 
элементов  из  облака  и  выяснение  химического  состава  облака,  а  также 
определение  его  первоначальной  массы.  Исследование  процесса  конденса- 
ции газа  в  пыль  позволило  бы  найти  относительное  содержание  газа  и  пыли 
и  определить  степень  непрозрачности  облака  и  его  температуру.  На  основа- 
нии этих  данных  уже  можно  было  бы  определить  скорость  уплощения  облака 
и  выяснить  роль  гравитационной  неустойчивости.  Кроме  того,  важно  было 
бы  исследовать  роль  турбулентных  движений  и  определить  среднюю  хаоти- 
ческую скорость  вещества  в  облаке,  роль  вязкости  и  скорость  передачи 
момента  количества  движения  из  центральных  областей  облака  во  внешние, 
роль  эффекта  Пойнтинга  —  Робертсона.  Необходимо  также   уточнение    ме- 
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ханизма  аккреции.  Многие  из  этих  вопросов  тесно  связаны  друг  с  другом, 
поэтому  для  строгого  решения  задачи  необходимо  комплексное  рассмотре- 
ние ряда  вопросов,  касающихся  эволюции  протопланетного  облака  и  обра- 
зования в  нем  планет. 

В  заключение  выражаю  благодарность  Б,  Ю.  Левину  -за  ряд  критиче- 
ских замечаний. 

Геофизический  институт  Поступила  в  редакцию 

Академии  наук  СССР  16  апреля  1954  г. 
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